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Abstract 

The Hodge equations for 1-forms are studied on Beltrami's pro- 
jective disc model for hyperbolic space. Ideal points lying beyond 
projective infinity arise naturally in both the geometric and analytic 
arguments. An existence theorem for weakly harmonic 1-fields, chang- 
ing type on the unit circle, is derived under Dirichlet conditions im- 
posed on the non-characteristic portion of the boundary. A similar 
system arises in the analysis of wave motion near a caustic. A class 
of elliptic-hyperbolic boundary-value problems is formulated for those 
equations as well. For both classes of boundary- value problems, an ar- 
bitrarily small lower-order perturbation of the equations is shown to 
yield solutions which are strong in the sense of Priedrichs. MSC2000: 
35M10, 58J32, 53A20, 78A05 



1 Introduction 

The projective disc was introduced by Beltrami^ in 1868. His construction 
was an early example of a Euclidean model for a non-Euclidean space, in 
this case, a space having curvature equal to —1. The projective disc has 
the striking property that even points infinitely distant from the origin are 
enclosed by the Euclidean unit circle centered at the origin of . This implies 
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the possibility of points in projective space which he beyond the curve at 
infinity. It is known that such ideal points arise naturally in the process of 
constructing normal and translated lines for chords of the projective disc. In 
this sense ideal points may be said to be intrinsic to the model, rather than 
only a theoretical possibihty allowed by the model. We call the union of the 
conventional projective disc and its ideal points the eoctended projective 
disc. 

Hua^ considered a second-order partial differential equation for scalar 
functions on the extended projective disc. He proved the existence of so- 
lutions to certain boundary-value problems of Tricomi type, in which data 
are given on characteristic curves, which represent trajectories of generalized 
wavefronts. Hua's work was extended to other problems of Tricomi type by 
Ji and Chen.^°'^^ The existence of a class of weak solutions to the Hodge 
equations for harmonic 1-fields on extended P^, with data prescribed only on 
the non-characteristic part of the boundary, was proven in Ref. 23. Locally, 
the Hodge equations reduce in the smooth scalar case to the equation studied 
by Hua. 

This communication provides a geometric and analytic context for such 
results (Sec. 1). In addition, we prove an existence theorem for weakly har- 
monic 1-fields which includes the results of Ref. 23 as a special case (Sec. 2.1), 
and consider a similar system that arises in optics (Sees. 3.1, 3.2). Boundary- 
value problems are formulated for both systems, in which the boundary con- 
tains points in both the elliptic and hyperbolic regions of the equations. 
These problems are shown in Sees. 2.2 and 3.3 to be an arbitrarily small, 
lower-order perturbation away from problems possessing a unique, strong 
solution. 

Because both scalar equations and systems are discussed, we distinguish 
a vector-valued solution by writing it in boldface. However, for typographic 
simplicity, coefficient matrices and operators are not written in boldface. 

1.1 A geometric classification of linear second-order 
operators 

The highest-order terms of any linear second-order partial differential equa- 
tion on a domain ^2 C can be written in the form 

Lu^a (x, y) u^^ + 2(5{x, y)u^y + 7(2;, y)uyy, 
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where {x,y) are coordinates on ft and a, (3, and 7 arc given functions. (A 
subscripted variable denotes partial differentiation in the direction of the 
variable.) 

If the discriminant 

A {x, y) =a-f- 0^ 

is positive, then the equation associated with the operator L is said to be of 
elliptic type. The simplest example is Laplace's equation, for which a — ^ — 1 

and /3 = 0. If the discriminant is negative, then the equation associated with 
the operator L is said to be of hyperbolic type. The simplest example is the 
normalized wave equation, for which q; = 1,7 = — l,/5 = 0; other forms are 
q; = — 1, 7 = 1, /? = 0, or q; = 7 = 0, /? = 1. If A = 0, then the equation 
associated with the operator L is said to be of parabolic type] examples are 
equations which model diffusion. If the discriminant is positive on part of Q 
and negative elsewhere on then the equation associated with the operator 
L is said to be of mixed elliptic-hyperbolic type. A simple example of an 
elliptic- hyperbolic equation is the Lavrent'ev-Bitsadze equation, for which 
a — sgn (y) , P = 0, and 7 = 1. 

If we take Q to be a smooth but curved surface, then we may not be 
able to cover Q hj a single system of Cartesian coordinates. However, we 
can always introduce Cartesian coordinates {x^,x'^) locally on any smooth 
surface, in the neighborhood of a point on the surface. In terms of such 
coordinates, the distance element ds on Q can be written in the form 

2 2 

ds^ = 9ij{x^, x'^)dx'^dx\ 

i=i j=i 

where gij is a symmetric 2x2 matrix, the metric tensor on Jl. (In the 
sequel we will understand repeated indices to have been summed from 1 to 
dim (n) without writing out the summation notation each time.) A natural 
differential operator on functions u defined on such a space is the Laplace- 

Beltrami operator 

where g'^^ is the inverse of the matrix g^j and g is its determinant. 

Laplace's equation can be associated to the Laplace-Beltrami operator 
on the Euclidean metric for which gij is the identity matrix. The wave 
equation for /3 = can be associated with the Laplace-Beltrami operator 
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on the 2-dimensional Minkowski metric gu — 1, (?22 = ~1) 5*12 = fi'21 = 0. 
The Lavrent'cv-Bitsadze equation can be associated to the Laplace-Beltrami 
operator on a metric which is Euchdean above the a;-axis and Minkowskian 
below the x-axis. 

In this classification, the type of a linear second-order equation is not a 
function of the associated linear operator at all; that operator is always the 
Laplace-Beltrami operator. Rather, the type of the equation is a feature of 
the metric tensor on an underlying surface. A Riemannian metric, in which 
the distance between distinct points of fl is always positive, corresponds 
to an elliptic equation, whereas a pseudoriemannian metric, for which the 
distance between distinct points may be zero, corresponds to a hyperbolic 
equation. The Laplace-Beltrami operator on a surface for which the metric 
is Riemannian on part of a surface and pseudoriemannian elsewhere will be 
of mixed elliptic-hyperbolic type. However, any sonic — or parabolic — curve 
on which the change of type occurs will necessarily represent a singularity of 
the metric tensor, as the determinant g will vanish along that curve. (The 
term sonic curve is borrowed from compressible fluid dynamics, in which the 
equations for the velocity fleld of a steady ideal flow change from elliptic to 
hyperbolic type at the speed of sound. The underlying pseudoriemannian 
metric in that case is called the flow metric.^) 

One definition of the signature of a metric is the sign of the diagonal 
entries of the metric tensor. Any change in the signature which results in a 
change in sign of the determinant g will change the Laplace-Beltrami operator 
on the metric from elliptic to hyperbolic type. The Laplace-Beltrami operator 
on surface metrics for which such a change occurs along a smooth curve will 
correspond to planar elliptic-hyperbolic operators in local coordinates. 

If we consider the distance element 

ds\ = a. (,x, y) dy^ — 2(5 (x, y) dxdy + 7 (x, y) dx'^, 

then null geodesies on the corresponding surface are solutions of the ordinary 
differential equation 

dsl = 0. 

The graphs of these solutions are called characteristic curves of the equation 
Lu — 0. Hyperbolic operators, which are associated with wave propagation, 
always have real-valued characteristics, or null geodesies. 

In determining the qualitative behavior of solutions to partial differential 
equations we often ignore lower-order terms, but this neglect is only justified 
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when considering purely second-order properties such as the nature of the 
sonic curve. The importance to this paper of lower-order terms is related 
to the fact that the Laplacc-Bcltrami equations on the extended projective 
disc are not of real principal type in the sense of Ref. 6; see Ref. 27 for an 
accessible discussion of scalar elliptic-hyperbolic operators of real principal 
type and their properties. 

1.2 The geometry and analysis of ideal points 

Here we review basic properties of Laplace-Beltrami equations on Beltrami's 
hyperbolic metric on the projective disc: 



(see , e.g., Ref. 32, Vol. I, Sec. 65 and Vol. II, Sec. 138, for a derivation). In 
this metric the unit circle is the absolute: the locus of points at infinity. 

The existence of points lying beyond the curve at infinity on the projective 
disc is natural from a geometric point of view. For example, choose a point 
p in the interior of the projective disc and draw a vertical line £^ through it. 
A hyperbolic line in the Beltrami metric is any open chord of the unit circle, 
so iy is a hyperbolic line plus two points at infinity and an ideal extension to 
points outside the unit circle. Denote by F{p) the family of hyperbolic lines 
created by rotating £y about p. Move p along the horizontal fine £h through 
p, and consider the affect of this motion on the family F[p). As p passes 
through the boundary of the unit circle k into the M^-complement of k, the 
family of hyperbolic rotations becomes a family of hyperbolic translations. 
For this reason, hyperbolic translations inside the unit disc can be interpreted 
as rotations about a point in lying beyond the unit disc. 

As another example, consider that the pole of a hyperbolic line i is the 
intersection of those two tangents to the unit circle which intersect i at the 
two points of its contact with the unit circle. (We call these the polar lines 
of £.) Thus any two hyperbolic lines £i and £2 are orthogonal if and only if 
the pole of £2 lies on the ideal extension of £1 and vice-versa. 

These and other geometric constructions on extended are described in 
more detail in Chapter 4 of Ref. 28. 

In order to sec that ideal points also arise naturally in analysis, consider 
the Laplace-Beltrami operator on the projective disc with Beltrami's metric. 



dx'^ + 2xydxdy + {1 — x'^) dy"^ 
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We have 



L [u] = (1 



u. 



'XX 



2xyuxy + (l - y^) u. 



yy 



+lower order terms] . 



The characteristics of the equation L[u] = satisfy the ordinary differen- 
tial equation 



where, as is conventional, we take 6 to be the angle between the radial vector 
and the positive Solutions of eq. (2) correspond geometrically to the 

family of tangent lines to the unit circle centered at the origin of M^. 

Thus the characteristic lines always include ideal points and wave prop- 
agation can only occur on regions composed of such points. 

The Laplace-Beltrami equations on extended come with a natural 
gauge theory in the following sense: The characteristic equation is obviously 
invariant under the projective group. So although the equations in the form 
in which we study them change type on the unit circle in M^, they are projec- 
tively equivalent to a system which changes type on any conic section. Note 
that whereas classical gauge theories are invariant under groups of Euclidean 
motions, which are inertial transformations, this kind of gauge invariance is 
with respect to a group of non-Euclidean motions, which are non-inertial. 
Also, the gauge theories which are familiar from particle physics act "up- 
stairs" on a fiber bundle of physical states. The transformation group under 
which the Laplace-Beltrami equations are invariant acts "downstairs" on the 
underlying metric, in the manner of the gauge group of general relativity. 
Indeed, analysis of wave motion on extended has certain similarities to 
the analysis of wave motion in the vicinity of a light cone (c./. Ref. 30). The 
time-like and space-like regions are inverted, and characteristic lines for the 
Laplace-Beltrami equation are analogous to the paths of photons. 



(l - y^) dx^ + 2xydxdy + (1 - x'^)dy^ = 0. 



(1) 



This equation has solutions 



a; cos 6* + ysinO = 1, 



(2) 
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2 Harmonic 1-fields on the extended projec- 
tive disc 



We can solve, instead of the Laplace- Beltrami equation, a system of two 
first-order equations of the form 

\9\'"^d,{g^^^\u^^Q, (3) 

\ (diUj - djUi) = 0, (4) 

where Ui = Ui{x^, x"^), i = 1,2. As in the second-order equation, g^j is a metric 
tensor on the underlying surface. Solutions u = {ui,U2) of this first-order 
system are (locally) harmonic 1-fields. Notice that if the scalar function 
ip{x^,x'^) satisfies v^^i = ui and ip^^ = U2, then (p satisfies the Laplace- 
Beltrami equations. But there are solutions if of the Laplace-Beltrami system 
for which the pair (<^a;i, (px'^) is not a harmonic 1-field. 

Consider a system of first-order equations on having the form 

Lu = f , (5) 

where 

L^{Lr,L2), f =(/i,/2), 
u = (mi {x, y) , U2 (x, y)) , (x, |/) G O CC R^. 
Let u satisfy (5) with 

(Lu)i = [(l - x^) ui] ^ - 2xyuiy + [{l - y'^) 1^2] ^ + hxui + k2yu2, (6) 

and 

(Lu)2 = (1 - y^) {uiy - U2x) + hxui kAyu2, (7) 

where VL is chosen so that 7^ 1 there. Here ki,k2,k3 and k^ are con- 
stants representing lower-order coefficients. In this section we consider three 
particular distributions of lower-order terms, studied in Ref. 23: 

Case 1: ki = /o^ = —2, k^ — k^ — 0. The domain of eqs. (5)-(7) in this 
case will be called Qi. 

Case 2: ki = —2, k2 = k^ = 0, k^ = 2. The domain of eqs. (5)-(7) in this 
case will be called Q2- 
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Case 3: ki = k2 = — — 0. The domain of cqs. (5)-(7) in this case 
wiU be called ^3. This case corresponds to eqs. (3), (4). 

The union of the domains Vti^ 1^2, and Vt^ will be called VL. 

A system of first-order equations can also be said to be of elliptic or 
hyperbolic type, and thus may change type along a singular curve. See, e.g., 
Ref. 5, Ch. III. 2. The higher-order terms of the preceding system can be 
written in the form A^Wx -\- A^Wy, where 





A^ 





(1 - y') 



(8) 



and 



—2xy I — y'^ 




(9) 



If 7^ 1, the characteristic equation 

\A' - AA^I = - (1 - y') [{1 - y') + 2xyX + (l - x')] 

possesses two real roots Ai, A2 on Q, precisely when x"^ + y^ > 1. Thus the 
system is elliptic in the intersection of Q with the open unit disc centered 

at (0, 0) and hyperbolic in the intersection of fl with the complement of the 
closure of this disc. The boundary of the unit disc, along which this change in 
type occurs, is the line at infinity on the projective disc and a line singularity 
of the tensor g^. 

Denote by a region of the plane for which part of the boundary dQ 

consists of a family of curves T composed of points satisfying eq. (1) and the 
remainder C = dQ\r of the boundary consists of points {x, y) which do not 
satisfy eq. (1). We seek solutions of eqs. (5)-(7) which satisfy the boundary 
condition 

dx dy „ 
ui— + U2-f- = 0, 10 

as as 

where s denotes arc length, on the non-characteristic part C of the domain 
boundary. Because the tangent vector T on C is given by 

^ dx . ^ dy . 
ds dx^^ 

a geometric interpretation of this boundary condition is that the dot product 
of the vector u = (1*1,1*2) and the tangent vector to C vanishes, i.e., u 
is normal to the boundary dfl on the boundary section C. We call these 
homogeneous Dirichlet conditions. 
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2.1 Weak solutions 



In Ref. 23, weak solutions to (5)- (7), (10) are shown to exist in certain 
weighted spaces on a class of domains. Here we extend that result to 
the case in which the domain is formed by the polar lines of a hyperbolic 
line I and a smooth curve C extending between the two polar lines of ^. 
The curve C must have the property that dy\c < when dVt is traversed in 
a counterclockwise direction. However, as long as this condition is met, C 
need not intersect the polar lines of £ at their points of tangency with the 
unit circle. Thus C may extend into both the elliptic and the hyperbolic 
regions. 

This domain is the analogue of the "ice-cream cone" -shaped domain as- 
sociated to the Tricomi equatiov?^ 

where in our case the curve C is the boundary of the ice-cream part and the 
polar lines, which are characteristics of eqs. (5)- (7), are the boundary of the 
cone part. The unit circle is the analogue of the x-axis, which is the sonic 
curve for the Tricomi equation. 

We initially consider the distribution of lower-order terms in case 1 of 
eqs. (6), (7). Let 6 lie in the interval [0,7r/4] and denote by Qi the region of 
the first and fourth quadrants bounded by the characteristic line 



the characteristic line 



Fi : X cos 9 + ysinO = 1, 



r2 : X cos ^ — y sin ^ = 1, 



and a smooth curve C. Let C intersect the hues Fi, r2 at two distinct 
points ci, C2, respectively. Assume that V(x, ?/) e Qi, l/-\/2 < X < V2 
and -l/v/2 < y < l/\/2, and that dy < on C A cusp is permitted for 
9 = 7r/4 at the points Ci, C2 = {1/^/2, ±l/\/2)- Otherwise, the boundary will 
have piecewise continuous tangent (so that Green's Theorem can be applied 
to it). Note that the domain considered in Sec. 3 of Ref. 23 is equivalent to 
this domain in the degenerate special case ^ = 0. 

Define U to be the vector space consisting of all pairs of measurable 
functions u = (^1,^2) for which the weighted norm 



u 




1/2 

(|2a;2 -l\uj+ |2|/2 - l| uj) dxdy 
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is finite. Notice that this expression vanishes at the intersection of £ with its 
polar fines at tfie value 9 = t^/A. Denote by W tfie linear space defined by 
pairs of functions w = (wi, W2) having continuous derivatives and satisfying: 



onr = riUr2; 

on C; and 



Widx + W2dy — 



(11) 
(12) 




Hi 



|2a;2-l| ^(L*w)i + |2y2- l| ^ {L*wf^ dxdy < 00. 



Here 



and 



(L*w), = [(1 - x^) wi]^ - 2xywiy + [(l - y^) W2]^ + 2xwi, 



(L*w)2 = (1 - y^) {wiy - W2x) + 2ywi. 



Define the Hilbert space H to consist of pairs of measurable functions h = 
{hi, /i2) for which the norm 




2x^ 



1-1 U2 



1 hi + \2y^-l\ 



^ hlj dxdy 



nl/2 



is finite. 

We say that u is a weak solution of the system (5)-(7), (10) in case 1 on 
ill if u e [/ and for every w e 1^, 



where 



-(w,f) = (rw,u), 



(w, f) = j J {wifi + ^2/2) dxdy. 



In case 2, we restrict the domain fl2 to lie in the fourth quadrant of 
the Cartesian plane. Define Fi, 12 C F to be characteristic fines which are 
tangent to the unit circle at distinct points in the fourth quadrant and which 
intersect at a point in the complement of the unit disc in R^. The curve C is 
defined analogously to the corresponding curve of fli. In particular, dy\c < 
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on Q2 when C is traversed in a counter-clockwise direction. Replace U by 
the space U' of all pairs u of measurable functions {ui,U2) for which the 
weighted norm 



u 



/ / (^xul + \y\ ufj dxdy 



1/2 



is finite. Replace W by the space W defined by pairs of continuously dif- 
ferentiable functions w = (tfi,W2) satisfying eq. (11) on F, eq. (12) on C, 
and 

I I [x'^ {L*w)\ + \y\~^ {L*wf^ dxdy < 00. 
J Jill 



In this case 



(L*w)i = [(1 - x^) wi]^ - 2xywiy + [(1 - y^) W2\y + 2xwi, 



and 



(L*W)2 = (1 - y'^) {Wly - W2x) - 2yw2. 

Finally, we replace H by the space H' of measurable functions h 
for which the norm 



(/il,/i2) 



ihir* = 



''"^1 + |y| ^ hl^ dxdy 



1/2 



is finite. 

Because is nonzero in case 2, the consistency condition (4) is violated 
and u cannot be the gradient of a scalar potential, even locally. Harmonic 
fields in which condition (4) is violated arise in various contexts — see Section 
4 of Ref. 25 for a nonlinear example — and correspond physically to stationary 
fields having sources. 

In case 3, we restrict the domain, ^3, to lie in the first quadrant. Define 
Fi, F2 C F to be characteristic lines which are tangent to the unit circle at 
distinct points in the first quadrant and which intersect at a point in the 
complement of the unit disc in R^. In this case we replace U and H hy L'^. 
We replace W by the space of pairs of functions {wi,W2) which satisfy 
(11) on F and (12) on C. Note that L is self-adjoint in case 3. In addition, 
we fix positive numbers S « 1/2 and s « 1/2 and require O3 to lie in the 
semi-infinite rectangle 
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1 



1 



<y < \/l - e. 




Weak solutions in cases 2 and 3 are defined exactly analogously to case 
1, with appropriate replacement of the domain and function spaces. 

Theorem 1. Let the lower-order terms in eqs. (6), (7) he distributed as 
in cases 1, 2, or 3, on the domains fli, fl2, or fl^, respectively. Then there 
exists a weak solution of the boundary-value problem (5)-(7), (10) for every 



Proof. The proof is an extension of the arguments in Ref. 23, so we will 
be brief. We derive a basic inequality, that there is a X e M"*" such that 
Vw e W, 



(with the norms appropriately adjusted in cases 2 and 3). We derive this 
inequality by choosing a scalar multiplier a, computing the inner product 
(L*w, aw) , and integrating by parts. Denoting the coefficients of wl off the 
boundary by a, those of wl by 7 and those of W1W2 by 2/3, we choose, in case 



K ||w||^ < ||L*w 



1, a = and obtain 



a — X 



and 



(3 = 



where 



2(3w\W2 > — 2x \xwi\ \yw2\ > — (• 
In case 2 we choose a — 1 and obtain 



'3 2, 2 2\ 

X W]^-\- xy W2) ■ 



a = 2a;, 7 = -2y, 



and /3 = 0. 

In case 3 we choose a — xy and obtain 



« = |(3^^-1),7 = |(1-Z/^), 
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and 



2/? = - (1 - X. 



The quadratic form — can be shown to be non-negative in case 3 by 
noticing that the argument in Sec. 6.2 of Ref. 23 does not use the restriction 
X < 1 and thus extends to OTir more general case. 

The remainder of the proof is essentially the same for all three cases. Ap- 
plying Green's Theorem to derivatives of products in (I/*w, aw) , we obtain 
a boundary integral / having the form 



Because wi vanishes identically on C, the boundary integral is nonnegative 
on C by the hypothesis on dy\c- On the characteristic curves, we no longer 
have the property that dx — 0, which we used in deriving the basic inequality 
of Ref. 23. However, 






Jr 



/^{(i-.>;.,+[2.,»;-(i-/M4. 



where we have used the fact that 



W2dy — —Widx 



on characteristic lines. In fact, we have 




( 



dy_ 
dx 



) 



-|- 2xyw\ — (1 — y'^)wiW2 dx 




by the same identity. Equation (1) implies that 




so we can write 
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C a 



2xy— + l-y 
ax 



WiW2dx+ 



2xywi{-W2^) - (1 - y'^)wiW2 



dx = 0. 



This establishes the basic inequahty. 

Proceeding as in Ref. 19, we use the basic inequality to apply the Riesz 
Representation Theorem and obtain an element h e for which 

-(w,f) = -(LV,h)*, 

where the product on the right is the inner product on H (or on H' or in 
cases 2 or 3, respectively). Writing hi and /i2 of h in terms of appropriate 
rescalings of ui and U2,^^ we obtain 



-(L*w,h)* = (L*w, u), 



which completes the proof. 



2.2 Strong solutions 

By a strong solution of the boundary-value problem (5), (10) we mean an 
element u G L^{Q) for which there exists a sequence u'' of continuously 
differentiable vectors satisfying the boundary condition (10), for which 



lim llu'^ — u| 



L2 



and 



0, 



0. 



lim IlLu'' 

For u = {ui{x,y),U2{x,y)), {x,y) e fl CC R^, define the operator L — 
(Li, L2) by the matrix equation 



Lu = A^n^ + A^Uy + Bu 



(13) 



for matrices A^, A^, and B. We say that L is symmetric-positive^ '^^'^^ if the 
matrices A^ and A^ are symmetric and the matrix 
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Q = 2B*-Al-Al 



is nonnegative. Here 

B*^^{B + B'), 

where for a matrix W = [wij], = [wji]. 

In cases for which L is not symmetric-positive, there may be a nonsingular 
matrix E such that EL is symmetric-positive. In that case we replace the 
equation 

Lu = f 

by the equation 

ELu = Ef 

and try to show that the operator EL is symmetric-positive. (The conversion 
of L into a symmetric-positive operator by the construction of a suitable 
multiplier E will not be used in this section, but will be used in Sec. 3.3.) 

Suppose that N{x,y), {x,y) G dfl, is a linear subspace of the vector 
space y, where u is regarded as a mapping u : flUdQ ^ V, and that N{x, y) 
depends smoothly on x and y. Define the matrix 

where n = (ni, 712) is the outward-pointing normal vector to dQ. The bound- 
ary condition that u lie in N is said to be admissible^^ if iV is a maximal 
subspace of V and if the quadratic form (u, /3u) is non-negative on dQ. 
A sufficient condition^ for admissibility is that there exist a decomposition 

P = P+ + P-, 

for which the direct sum of the null spaces for (3j^ and spans the restriction 
of V to the boundary, the intersection of the ranges of /9+ and (3- have only 
the vector u = in common, and the matrix = /?+ — /3_ satisfies 

^ 2 - 

In this case the boundary condition 

p_u = Qondn 
15 



is admissible for the boundary- value problem 

Lu = f inVL. 
Moreover, the boundary condition 

/5^w = QondVL 
is admissible for the adjoint problem 

L*w — YiinQ.. 

These two problems possess unique, strong solutions whenever the differ- 
ential operators are symmetric-positive and the boundary conditions are 
admissible/'^^ 

In this section we give sufficient conditions for the existence of certain 
strong solutions arising from an arbitrarily small lower-order perturbation of 
the Laplace-Beltrami equations on extended P^. We do so by showing that 
the differential operator L given by (5)- (7) with fci = A;2 = ^3 = ^4 = is 
arbitrarily close to a symmetric-positive operator and by stating an admis- 
sible boundary condition. The existence of strong solutions to a different 
perturbation on an explicit domain will be shown in Sec. 3.3. 

If the matrices and of eq. (13) are given by eqs. (8) and (9) and 
the matrix B is given by 

/ -2x -2y 
\ 

then the quantity Q is zero. Thus we replace the matrix S by a matrix 
which differs from B by an arbitrarily small perturbation and takes the form 

/ -2x + ei -2y + e2 \ .^.x 

where ei > 0, £4 > 0, £2 + (1 — y^) £3 > 0, and 

[e2 + {l-y')esY<Hl-y') e,e,. 

If we choose the domain of L in such a way that y"^ < 1 there, then this re- 
placement converts Q into a positive-definite matrix and L into a symmetric- 
positive operator. 
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Denote by fl^ a domain having boundary — E U F such that 
y"^ < 1 on Q4. Let the components of the normal vector n on 8^4 be given 
by (ni,n2). Assume that rii and n2 never vanish at the same point of 8^4. 
We place conditions on ni, n2, and dfl4 sufficient to guarantee admissibility 
of the boundary condition 

1*1^2 — U2ni — (15) 

on F, with no condition given on E. 

Let ni > 0, n2 < on F and ni < 0, n2 > on E. Defining the adjoint 
space as in Sec. 2.1, for w e V* we take w = (0, W2) on F and 

Win2 — W2ni — 

on E. Define 

a = [- (1 - y^) n2/ni + 2xy - (l - x'^) rii/n-^ n2. 

Assume that a — QonE, and that a < on F. 

Theorem 2. The boundary-value problem 

LvL = A^n^ + A^xxy + B^n = f 

for {x,y) e with , A'^, and given by eqs. (8), (9), and (I4) respec- 
tively and with condition (15) imposed on the curve F of 80,4, possesses a 
unique, strong solution u{x,y) for every f e L^(Jl4). 

Proof. Because the matrix B^ has been constructed in such a way that L 
is symmetric-positive, it remains only to show that the boundary condition 
(15) is admissible on ^4. 

We have 

f -a-{l-y'^)nl/ni {l-y'^)n2 \ 
{l-y')n2 -{l-y')mj- 

Note that the apparent singularities in /3 at ni = and in a at ^2 = are 
removable. 

On F, choose 
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and 



o (a 2\ [ -^2/^1 

/3_ = (1 - 7/ ) 712 ( ^ 

On choose 

a 2\ ( -n2/ni 

= (1 - y j ^2 ( ^ 

and 

If u G Vji?, then (15) imphes that /3_u = 
that w*/3+ = for w G V*p. If u G V\e^ then /9_u = for all values of u and 
w*/3_|_ = by the properties of and a. So the direct sum of the null spaces 
of /3_ and /9+ spans V on dQ,^. Moreover, the hypotheses guarantee that the 
ranges of /3_ and /3+ have only the zero vector in their intersection. Finally, 

(5- = n* >Q 

on both E and F. 

This completes the proof of Theorem 2. 

3 An analogous problem from optics 

Geometrical optics is a zero- wavelength approximation to classical wave me- 
chanics in which the governing differential equations are replaced by the Eu- 
clidean geometry of rays. The limitations of the geometrical optics approxi- 
mation arc apparent in the neighborhood of caustics, which are envelopes of a 
family of rays. It is not simply that geometrical optics predicts infinite inten- 
sity in such regions, whereas diffractive effects reduce the predicted intensity 
to a finite number. Even in applications for which the agreement between 
the predictions of geometrical optics and experiment is generally good, the 
former may predict singularities, e.g., cusps, which are entirely smoothed 
out by diffraction. A dramatic example of this for the case of water waves 
is illustrated in Figures 5.6.1 and 5.6.2 of Ref. 29. This is, of course, far 



1 

-ni/n2 



1 

-ni/n2 



. The properties of V* imply 
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from the only drawback of the geometrical optics approximation. See, for 
example, the discussion of the rainbow caustic in Sec. 6.3 of Ref. 22. 

The accuracy of the geometrical optics approximation can be improved by 
considering waves of arbitrarily high frequency obtained by uniform asymp- 
totic approximation of solutions to the Helmholtz equation (Sec. 3.1). While 
the older of these approximations also fail at caustics, an asymptotic formula 
introduced independently by Kravtsov^^ and Ludwig^^ retains its meaning 
even in the neighborhood of a caustic; see Ref. 13 for a review. 

Recently, Magnanini and Talenti studied a nonlinear elliptic-hyperbolic 
equation, implied by the Ludwig-Kravtsov approximation, having the form^^ 

(iVt'l'' - Vy) Vxx + '2-VxVyVxy + (iVvl"^ - v^) Vyy = 0, (16) 

where v — v{x,y), {x, y) e R^. Those authors were able to show the existence 
of weak solutions to the full Dirichlet problem for the linear eUiptic-hyperbolic 
equation 

which is related to eq. (16) by the Legendre transformation 

VL{p,q) = xp + yq-v{x,y). (18) 

Magnanini and Talenti 's result is remarkable in that it is difficult to formu- 
late a full Dirichlet problem which is well-posed for a given elliptic-hyperbolic 
equation, even in the weak sense; by full we mean that data arc prescribed on 
the entire boundary. Morawetz's proof of the existence of weak solutions to 
the full Dirichlet problem for the Tricomi equation, the most intensively stud- 
ied elliptic- hyperbolic equation, required a delicate argument. The full 
Dirichlet problems for other important elliptic-hyperbolic equations remain 
unknown. For example, the full Dirichlet problem has not been correctly 
formulated even for weak solutions to a scalar elliptic-hyperbolic equation 
associated to electromagnetic wave propagation in cold plasma, although a 
well-posed Dirichlet problem for weak solutions has been formulated for data 
prescribed only on part of the boundary. (In fact, Magnanini and Talenti do 
more than prove the existence of a weak solution: they also show uniqueness 
and internal regularity modulo a point, and give an explicit representation 
of the solution in terms of special functions.) 
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The existence of a well-posed Dirichlet problem is important because 
physical reasoning often suggests that the full Dirichlet problem is the correct 
problem even in the case of equations for which mathematical reasoning 
suggests otherwise. 

Two questions suggested by Magnanini and Talenti's paper are: 
i) The transformation (18) itself fails at caustics (which are not generally 
identical to the caustics of the physical model). One would like to charac- 
terize regions at which this linearization method fails and the nature of the 
singularities that arise in such regions. See, for example, Proposition 2 of 
Ref. 26. 

a) The result proven in Ref. 18 requires the domain boundary to lie en- 
tirely within the elliptic region of the equation. It is an important quality of 
eq. (17) that the elliptic region surrounds the hyperbolic region, a property 
not shared by other elliptic-hyperbolic equations. Thus there is some math- 
ematical interest in asking whether solutions of (17) exist with boundary 
points lying in both the elliptic and hyperbolic regions, a situation in which 
this special condition is no longer applicable. We consider this question in 
Sec. 3.3. 

Equation (16) is a special case of the system 

(p' + g'] P. + 2pgp, + [(p' + p'] Qy = 0, (19) 

(20) 



Vy - Qx 



0. 



This system is equivalent to eq. (16) if there is a continuously differen- 
tiable scalar function v (x, y) for which = p and Vy = q. (Such a function 
always exists locally, by eq. (20).) 

Consider any two-dimensional quasilinear system of two equations having 
the form 



an ai2 

021 ^22 



d_ 
dx 



bn bi2 

b21 b22 



d_ 
dy 



(21) 



where the entries of the coefficient matrices depend only on p and q. Then 
the coordinate transformation {x, y) — > (p, q) takes eq. (21) into the linear 
form 



bi2 

b22 



-ai2 

-0.22 



d_ 
dp 



+ 



-bu an 
-621 O21 



d_ 
dq 
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provided the Jacobian of the transformation 

J d {x, y) dx dy dy dx 
d (p, q) dp dq dp dq 

is nonzero. This special case of the Legendre transformation is called a hodo- 
graph map, and the space having coordinates (p, q) is called the hodograph 
plane; see, e.g., Sec. V.2.2 of Rcf. 5. 

The coordinate systems {p,q) and {x,y) are related by eq. (18), where 

, , fdV dV 



dp dq 



and 



/ dv dv^ 
\dx^ dy ^ 

Applying a hodograph transformation to eqs. (19), (20) yields the system 



/ 2 I 2\2 2 

[p +q) -p 



Xp - 2pqxq + 



{p' + q^f-q^ 7/, = 0, 



(22) 

x,-yp = 0. (23) 

This system is equivalent to eq. (17) if there is a continuously differentiable 
scalar function V {x, y) for which Vp = x and Vq = y. (Again, this can always 
be arranged locally.) 

As in Sec. 2, we write the second-order terms of eqs. (22), (23) in the 
form A^Wx + A^Uy, where u = u(x, y) and in this case 



A^ 



{x' + y'f 




and 



X 







-2xy (x'^ + y^) - y^ 
1 



The characteristic equation 

_XA'\^-{ [{x' + y'Y - y'] + 2xyX + [{x' + y')' - x'] } 

possesses two real roots Ai, A2 precisely when x^ + y^ > {x"^ + |/^)^ , that is, 
when x^ +y^ < 1. Thus the system is hyperbolic at points lying inside the 
open unit disc centered at {x, y) — (0, 0) and elliptic outside the closure of 
this disc. The circle x'^ + y'^ = 1, along which the change in type occurs, is 
the parabolic region of the system. 
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3.1 Uniform asymptotic approximations 

Substitution of the simplest formula for an oscillatory wave into the wave 
equation results in the Helmholtz equation 

AU (x) + k'^i^'^U (x) = 0, (24) 

where we take x to be a vector in R^, and where A; and u are physical con- 
stants. In the standard application, u is the refractive index of the medium 
and k is inversely proportional to wavelength. In the region of visible light, 
the wavelength is sufficiently small that k dominates over all other mathe- 
matically relevant parameters, an undesirable property known as stiffness. 

For this reason, short-wave solutions of (24) are usually approximated 
by uniform asymptotic expansions^"^'^^ which satisfy (24) to arbitrarily high 
order in These approximations are valid in regions which contain smooth 
and convex caustics such as a circular caustic. The size of the region of 
validity is independent of k. Take u = 1 and approximate the solution to 
(24) by an expansion having the form 

U approxip^ 1 

\^Z{k-l'u) (^f:W^,(r)-(z^)-^j +^^'(^^/^«) {^X,{v).m-^^ 

X exp [ikv (x, y)] , 

where u (x, y) , v (,t, y) , Wj (r) , and Xj (r) are functions which do not depend 
on k and which are to be determined with the solution; the function Z{t) is 
a solution of the Airy equation 

Z" it) 

with initial conditions 

Z(0) 

and 

Z'{0) 

where F ( ) is the gamma function. 



- tZ {t) = 0, 

" F (2/3) 
F(l/3)' 
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This model implies the following system of equations for u and v : 

UxVx + UyVy = 0. 

In Ref. 18 three possible solutions of this system are enumerated: 

M = 0, \Vv\^ = 1; 

\Wu\ = 0, \Wvf = 1; 

the third possibility is that eq. (16) is satisfied. 

Obviously, the third alternative is the most interesting, and this case 
is studied in Ref. 18. This case is linearized to eq. (17) by a hodograph 
transformation. 

3.2 A first-order system 

Thus we are led to a system resembling eqs. (5)-(7): 

Lu = g, (25) 

where 

L = (Li,L2) , g = (91,92) , 
u = (mi (x, y) , U2 {x, y)) , (x, y) eVLdC 
(Lu)i = [/ {x, y) - x^] - '2xyuiy + [f (x, y) - y'^] U2y (26) 

and 

(Lu)2 = [/ (x, y) - y^] {uiy - U2x) , (27) 

for 

f{x,y)^{x' + y')\ (28) 
The domain is chosen so that 

/(a;,2/)-2/VO, 

under which system (25)- (28) becomes an inhomogeneous generalization of 
eqs. (22), (23). If in particular, 91 = 92 = 0, ui = Vx, and U2 = Vy, where 
V {x,y) is a scalar function, then eqs. (25)-(28) reduce to eq. (17). 
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As in the preceding sections, the second-order terms of eqs. (25)- (28) can 
be written in the form A^u^ + A^Uy, where 



A' 



f{x,y)-x^ 

-{f{x,y)-y') 



and 



-2xy f{x,y)-y'^ 
f{x,y)-y^ 

We find that the system is hyperbohc in the intersection of Q with the open 
unit disc centered at (0, 0) and elhptic in the intersection of Q with the 
complement of the closure of this disc. 



3.3 Strong solutions in an annulus 

Writing eq. (17) in polar coordinates {r,9), r > 0, < 9 < 2tt, we obtain^^ 

{r^-l)Vrr + rVr + Vee^O. (29) 

Letting ui — Vr and U2 — Vo transforms eq. (29) into a first-order system of 
the form 



Lu = A^Ur + A^ue + Bu = i, 
with u = (ui(r, 9),U2{r, 9)) J ^ (0, 0), 



A'^ 





-1 



,A' 



1 

1 ; ' 



(30) 



(31) 



and 



B = 



r 




As in Sec. 2.2, the matrices are symmetric and we find that Q = 2B* — 
Al — Al is exactly zero, suggesting that an arbitrarily small perturbation of 
the matrix B will result in a symmetric-positive operator. However, we find 
that we can retain the consistency condition Uig — U2r = if we employ a 
multiplier E as described in Sec. 2.2. Thus we define 



E 



a c(l — r^) 
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where a = a{r, 6) and c = c(r, 9) are continuously differentiable functions to 
be chosen. We replace B by the matrix 

where £i,£2 are arbitrarily small, positive constants. 
Replacing eq. (30) by the system 

EL = EA^Ur + EA^ug + EB,u = Ei, (33) 

with A^, A'^, and given by eqs. (31) and (32), we find that EL is a 
symmetric-positive operator provided we choose < Eq < r < R < oo, c a 
positive constant, and 

. (V2-£i)c 

where M is a constant such that M » c. 

We will solve eqs. (33) in the annulus Vt^ given by £o < r < V2. (The 
solutions can be patched into an elliptic boundary-value problem on the 
annulus \/2 < r < R.) Data will be prescribed on the outer boundary only. 
Annular domains are natural when numerical methods are used to study an 
equation, such as eq. (17), which is known to be singular at the origin, with 
the singular point excluded. The problem is also of some historical interest. 
An equation differing from (17) only in its lower-order terms was one of the 
first elliptic-hyperbolic equations to be studied, more than 75 years ago, by 
Bateman (Sec. 9 of Ref. 1). That equation equation arose from the solution of 
Laplace's equation in toroidal coordinates.^ At the time, Bateman raised the 
question of the existence and uniqueness of solutions in an annular region 
containing the unit circle, in which the outer boundary lies in the elliptic 
region and the inner boundary lies in the hyperbolic region of the equation. 

Although the system that we consider is a small perturbation of the one 
studied in Ref. 18, we note that the original equation is itself an approxima- 
tion, as described in Sec. 3.1. 

Theorem 3. Equations (33) with boundary conditions 

r{e)ui + a{e)u2 = 0, a^e) > r^e) (34) 
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imposed on the outer boundary r — -\/2, possess a strong solution on the 
annulus for every /e L'^{fl5). 



Proof. Although the equations are different, the argument is similar to 
the proof by Torre^° of the corresponding assertion for the hehcally reduced 
wave equation. 

The matrices E and have been constructed in such a way that the op- 
erator EL is manifestly symmetric-positive (for large M) , and the proof again 
reduces to a demonstration that the boundary conditions are admissible. At 
the outer boundary, choose Uouter — dr. Then 



outer 



Choose 



Then 



g-2 _|_ ^2 



(t2 + t2 



outer ~ 



arc + T a a c + ara 
—ara + r'^c —a'^a + arc 



P. 



'outer+ 



—arc + a a 
ara + a^c 



ara 
- arc 



Notice that Pouter+ + Pouter- = Pouter and that Pouter-'^ = 0, as (34) implies 
that U2 — — (t/(7)wi on the circle r = V^. Moreover, 



1 / (a2 - r2) a - 2aTC (r^ -a^)c- 2aTa 



a^ + r^ y {a^ -t^)c + 2ara (cr^ - r^) a - 2arc 
implying that 

* 

1 / (ct2 - t2) a - 2arc 



0-2 + t2 V (o-2 - t2) a - 2(7TC 

But this matrix is non-negative, given that > t^, provided that we choose 
M sufficiently large. 
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On the inner boundary we choose 



Then 



Choose 



Pi. 



a c 








Then Pinner+ = Pinner and /i* > foi M Sufficiently large. Moreover, 
Pinner-'^ = ou the circlc v = Eq for any vector u. 
This completes the proof of Theorem 3. 



4 A remark on terminology and notation 

Hodge* originally considered a p-form cu to be harmonic if it satisfies the 
first-order equations 

duj = Suj = 0, (35) 

where d : — >■ A^"*"^ is the exterior derivative and 6 : A^"*"^ — >■ A^ is the 
adjoint of d. If the underlying space is and a; is a 1-form given by 

LU = pdx + qdy, 

where p and q arc continuously differentiable functions, then the Hodge equa- 
tions (35) reduce to the Cauchy-Riemann equations for p and —q. However, 
although d is independent of the underlying metric, its adjoint 6 has a differ- 
ent local form for different metrics. Thus for a surface having metric tensor 
Qij, the Hodge equations for 1-forms are equivalent to the system (3), (4). A 
discussion of exterior forms and their properties is given in, e.g., Ref. 21. 

The standard definition of a harmonic form is given in terms of a second- 
order operator: it is a solution of the form- valued Laplace-Beltrami equations 

{d5 + Sd)u;^ 0. 

If the domain has zero boundary (either no boundary or the prescribed value 
a; = on the boundary), then the definitions in terms of first- and second- 
order operators are equivalent. Otherwise, one distinguishes them by calling 
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a form that satisfies eqs. (35) a harmonic field. In words, the Hodge equations 
assert that a harmonic field uj is both closed {duj = 0) and co- closed {6uj = 0) 
under the exterior derivative d. Obviously, every harmonic field is a harmonic 
form, but the converse is false. 

Notice that in eqs. (6) and (7), Li ^ 5 and L2 ^ d. For example, L-2 
includes a factor of 1 — whereas d does not, and 5 includes determinants 
of the metric tensor, whereas Li does not. In addition, cases 1 and 2 of 
(6), (7) include additional lower-order terms. Thus for example 5 and d are 
self-adjoint, whereas Li and L2 are not unless ki = k2 = = = 0. 

Acknowledgment. I am grateful to an anonymous referee for helpful 
criticism of an earlier draft of this paper. 

5 References 

^Bateman, H., "Notes on a differential equation which occurs in the two- 
dimensional motion of a compressible fluid and the associated variational 
problems," Proc. R. Soc. London Ser. A, 125, 598-618 (1929). 

^Bateman, H. Partial Differential Equations., Dover, New York, 1944. 

^Beltrami, E., "Saggio di interpretazione della geometria non-euclidea," Gior- 
nale di Matematiche 6, 284-312, 1868. 

^Bers, L., Mathematical Aspects of Subsonic and Transonic Gas Dynamics, 
Wiley, New York, 1958. 

^Courant, R. and D. Hilbert, Methods of Mathematical Physics, Vol. 2, 
Wiley-Interscience, New York, 1962. 

^Dencker, N., "On the propagation of polarization sets for systems of real 
principal type," J. Functional Analysis 46, 351-372 (1982). 

^Priedrichs, K. O., "Symmetric positive linear differential equations," Com- 
mun. Pure Appl. Math. 11 333-418 (1958). 

^Hodge, W. V. D., "A Dirichlet problem for harmonic functionals with appli- 
cations to analytic varieties," Proc. London Math. Soc. 36, 257-303 (1934). 



28 



^Hua, L. K., "Geometrical theory of partial differential equations," in: Pro- 
ceedings of the 1980 Beijing Symposium on Differential Geometry and Differ- 
ential Equations (S. S. Chern and Wu Wen-tsiin, eds.), Gordon and Breach, 
New York, 1982, pp. 627-654. 

^°Ji, X-H., and D-Q.Chen, "Tricomi's problems of non- homogeneous equation 
of mixed type in real projective plane," in: Proceedings of the 1980 Beijing 
Symposium on Differential Geometry and Differential Equations (S. S. Chern 
and Wu Wen-tsiin, eds.), Gordon and Breach, New York, 1982, pp. 1257- 
1271. 

^^Ji, X-H., and D-Q. Chen, "The Tricomi's problem of the nonhomogeneous 
equation of mixed type in the real projective plane," in: Mixed Type Equa- 
tions (J. M. Rassias, ed.), Teubner, Leipzig, 1986, pp. 280-300. 

^^Kravtsov, Yu. A., "A modification of the geometrical optics method" [in 
Russian], Radiofizika 7 664-673 (1964). 

^^Kravtsov, Yu. A. and Yu. I. Orlov, Gaustics, Gatastrophes, and Wave 
Fields, Springer- Verlag, New York, 1999. 

^"^Ladyzhenskaya, O. A., and N.N. Ural'tseva, Linear and Quasilinear Elliptic 
Equations, Academic Press, New York, 1968. 

^^Lax, P. D., and R. S. Phillips, "Local boundary conditions for dissipative 
symmetric linear differential operators," Gommun. Pure Appl. Math. 13, 
427-455 (1960). 

^^Lin, C. S., "The local isometric embedding in of two-dimensional Rie- 
mannian manifolds with Gaussian curvature changing sign cleanly," Gom- 
mun. Pure Appl. Math. 39 867-887 (1986). 

^^Ludwig, D., "Uniform asymptotic expansions at a caustic," Gommun. Pure 
Appl. Math. 19, No. 2, 215-250 (1966). 

^^Magnanini, R., and G. Talenti, "Approaching a partial differential equation 
of mixed elliptic-hyperbolic type," in Ill-posed and Inverse Problems, S. I. 
Kabanikhin and V. G. Romanov, eds., VSP 2002, pp. 263-276. 

^^Morawetz, C. S., "A weak solution for a system of equations of elliptic- 
hyperbolic type," Gommun. Pure Appl. Math. 11, 315-331 (1958). 



29 



^°Morawetz, C. S., "The Dirichlet problem for the Tricomi equation," Com- 
mun. Pure Appl. Math. 23, 587-601 (1970). 

^■"^Morita, S., Geometry of Differential Forms, American Mathematical Soci- 
ety, Providence, 2001. 

^^Nye, J. F., Natural Focusing and the Fine Structure of Light, Institute of 
Physics Pubhshing, Bristol, 1999. 

^^Otway, T. H., "Hodge equations with change of type," Annali Mat. Pura 
ed Applicata 181, 437-452 (2002). 

^^Otway, T. H., "A boundary- value problem for cold plasma dynamics," J. 
Appl. Math. 3, 17-33 (2003). 

^^Otway, T. H., "Maps and fields with compressible density," Rendiconti 
Sem. Mat. Univ. Padova 111, 133-159 (2004). 

^^Otway, T. H., "Geometric analysis near and across a sonic curve," in: New 
Developments in Mathematical Physics Research (C. V. Benton, ed.). Nova 
Science Pubhshers, New York, 2004, pp. 27-54. 

^^Payne, K. R., "Interior regularity of the Dirichlet problem for the Tricomi 
equation," J. Math. Anal. Appl. 199, 271-292 (1996). 

^^Stillwell, J., Geometry of Surfaces, Springer, New York, 1992. 

^^Stoker, J. J., Water Waves, Interscience, New York, 1957. 

^°Torre, C. G., "The helically-reduced wave equation as asymmetric-positive 
system," J. Math. Phys. 44, 6223-6232 (2003). 

Tricomi, F., "SuUe equazioni lineari alle derivate parziali di secondo ordine, 
di tipo misto," Rendiconti Atti dell' Accademia Nazionale dei Lincei, Ser. 5, 
14, 134-247 (1923). 

^^Veblen, O, and J. W. Young, Projective Geometry, Ginn and Co., Boston, 
1918. 



30 



